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Abstract 

Electric fields in non-local media with power-law spatial dispersion are discussed. 
Equations involving a fractional Laplacian in the Riesz form that describe the electric 
fields in such non-local media are studied. The generalizations of Coulomb’s law and 
Debye’s screening for power-law non-local media are characterized. We consider simple 
models with anomalous behavior of plasma-like media with power-law spatial dispersions. 

The suggested fractional differential models for these plasma-like media are discussed to 
describe non-local properties of power-law type. 
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1 Introduction 


Fractional calculus is dedicated to study the integrals and derivatives of any arbitrary real 
(or complex) order. It has a long history from 1695 [H [2]. The hrst book dedicated specih- 
cally to study the theory of fractional integrals and derivatives and their applications, is the 
book by Oldham and Spanier |3] published in 1974. There exists a remarkably comprehensive 
encyclopedic-type monograph by Samko, Kilbas and Marichev [1]. Many other publications 
including different fractional models had appeared from 1980. For example, see 0 - 1221 . Frac¬ 
tional calculus and the theory of integro-differential equations of non-integer orders are powerful 
tools to describe the dynamics of anomalous systems and processes with power-law non-locality, 
long-range memory and/or fractal properties. 

Spatial dispersion is called the dependence of the tensor of the absolute permittivity of the 
medium on the wave vector [231 [241 [25]. This dependence leads to a number of phenomena, 
such as the rotation of the plane of polarization, anisotropy of cubic crystals and other [26l [23 
[251 [22115IJ1 [5n [521 [551 [541 [551 [551 [57] . The spatial dispersion is caused by non-local connection 
between the electric induction D and the electric held E. Vector D at any point r of the 
medium is not uniquely dehned by the values of E at this point. It also depends on the values 
of E at neighboring points r', located near the point r. 

Plasma-like medium is medium in which the presence of free charge carriers, creating as 
they move in the medium, electric and magnetic helds, which signihcantly distorts the external 
held and the ehect on the motion of the charges themselves [251 EH ESj. The term ’’plasma¬ 
like media” refers to media with high spatial dispersion. These media are ionized gas, metals 
and semiconductors, molecular crystals and colloidal electrolytes. The term ” plasma-like me¬ 
dia” was introduced in 1961, by Viktor P. Silin and Henri A. Rukhadze in the book ’’The 
electromagnetic properties of the plasma and plasma-like media” [25] . 

In Section 2 the basic concepts and well-known equations of electrodynamics of continuous 
media are considered to hx the notation. In Section 3, we consider power-law type generaliza¬ 
tions of Debye’s permittivity and generalizations of the correspondent equations for electrostatic 
potential by involving the fractional generalization of the Laplacian. The simplest power-law 
forms of the longitudinal permittivity and correspondent equations for the electrostatic poten¬ 
tial are suggested. The power-law type deformation of Debye’s screening and Coulomb’s law 
are discussed. These suggested simple models allows us to demonstrate new possible types of 
an anomalous behavior of media with fractional power-law type of non-locality. In Section 4 
the description of weak spatial dispersions of power-law type in the plasma-like media is dis¬ 
cussed. The fractional generalizations of the Taylor series are used for this description. The 
correspondent power-law deformation of Debye’s screening and Coulomb’s law are considered. 
A short conclusion is given in Section 5. In Appendix 1, we suggest a short introduction to 
the Riesz fractional derivatives and integrals. In Appendix 2, the fractional Taylor formulas of 
different types are described. 

In this section we review some basic concepts and well-known equations of electrodynamics 
of continuous media, to hx the notation. For details see [2511211 [25] . 
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The behavior of electric helds (E, D), magnetic helds (B, H), charge density p, and current 
density j is described by the well-known Maxwell’s equations 


divD(t,r) =p(t,r), 


curlE(t, r) = 


9B(t, r) 


dt ’ 
divB(t,r) = 0, 

5D(t, r) 


curl H(t, r) = j(t, r) -|- 


dt 


( 1 ) 

( 2 ) 

(3) 

(4) 


The densities p(t,r) and j(t,r) describe an external source of held. We assume that the 
external sources of electromagnetic held are given. The vector E(f, r) is the electric held 
strength and the vector D(t,r) is the electric displacement held. In free space, the electric 
displacement held is equivalent to hux density. The held B(t, r) is the magnetic induction and 
the vector H(t, r) is the magnetic held strength. 

In the case of the linear electrodynamics the constitutive equations (material equations) are 
linear relations. For electromagnetic helds which are changed slowly in the space-time, we have 
the constitutive equations (material equations) in the well-known form 


Di{t,r) = 


(5) 

( 6 ) 


where eij and pij are second-rank tensors. For helds varying in space rapidly, we should consider 
the inhuence of the held at remote points r' on the electromagnetic properties of the medium 
at a given point r. The held at a given point r of the medium will be determined not only the 
value of the held at this point, but the held in the areas of environment, where the inhuence of 
the held is transferred. For example, it can be caused by the transport processes in the medium. 
Therefore, we should use non-local space relations instead of equations (El), ([6]). These non-local 
relations take into account space dispersion. For linear electrodynamics we have the following 
relation between electric helds E and D given by 


A(t,r)= / eij{r,r') Ej{t,r') dr\ 


(7) 


and 


Bi{t,r) = / Hj{t,r')dr'. 


( 8 ) 


If the medium is not limited in space and homogeneous, then the kernel of the integral 
operator is a function of the position diherence r — r'. 


Di{t, r) = / iij{r - r') Ej{t, r') dr'. 


(9) 
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In this case we can use the Fourier transform. The direct and inverse Fourier transforms T 
and for suitable functions, are given by 

(.F/)(k)=7-[/(r)](k)= f e-‘'‘‘'>/(r)dV, (10) 

dR3 

^(r) = (-F"^^)(r) = J^~i[^(k)](r) = f g{k) d^k. (11) 

(zvr)'’ J^3 

We does not use the hat for D(f, r) and E(f, r) to have usual notation. From the context, it will 
be easy to understand if the field is considered in the space-time of its the Fourier transforms. 

Then electric field will be represented as a set of plane monochromatic waves, for which 
space-time dependence are defined by the function exp{ikr — icot}. Therefore, relation dH]) has 
the form 

Di{u, k) = eij{k)Ej{u, k). (12) 

The function £p (k) is called the tensor of the absolute permittivity of the material: 

e:ij(k) = [ e“**"“'dij(r) dr'. (13) 

dR3 

Even for an isotropic linear medium, the dependence of the tensor e:ij(k) of the wave vector 
k preserves tensor form [23112H 133] • In this case we have 

eij{k) = £±(|k|) + |^£||(|k|), (14) 

where £x(|k|) - the transverse permittivity, and £||(|k|) - the longitudinal permittivity. 

The Maxwell’s equations for the electromagnetic fields have the well-known form [2311241125] 

i(k,E(u;,k))£(k) = p(n;,k), (15) 

[k, E(a;, k)] = a; B(c(;, k), (16) 

(k,B(n;,k)) = 0, (17) 

[k,B(a;,k)] = -ia;£(k) E(a;, k) +j(n;,k). (18) 

In these equations we have neglected the frequency dispersion. This can be done when the 
inhomogeneous field can be approximately regarded as static. 

In the case of a static external field sources in the environment can create a inhomogeneous 
electric field E(f, r) = E(r). The electric field in the medium is given by 

E(r) = — grad <h(r), (19) 
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where *h(r) is a scalar potential of electric field. Relation flT^ . applying the Fonrier transform, 
can be written by 

E(k) = -?k<Fk. (20) 

Therefore, snbstitnting (|20|) into flT^ . we obtain 

|kp£||(|k|)«Fk = Pk, (21) 


where pk = p(0,k). Note that equation fl2lD does not depend of the transverse permittivity 
£x(|k|). 

When the held source in the medium is the resting point charge, then the charge density is 
described by delta-distribution 

p(r)=QiW(r). (22) 

Therefore the electrostatic potential of the point charge in the isotropic medium, according to 
the equation fl2T]) . has the form 


<h(r) 


Q 

(27r)3 


^+ik(r 


|kPe||(|k|) 


d\ 


(23) 


where $(r) is the electric potential created by a point charge Q at a distance |r| from the 
charge. 

Let us note the well-known case [2211211122] is such that £||(|k|) = Eq, where the constant Eq 
is the vacuum permittivity (^q ~ 8.854 • m~^). Substituting £||(|k|) = Eq into (I2T]) . we 

obtain 

|kp<Fk = -pk. (24) 

^0 

The inverse Fourier transform of fl2T|) gives 

A$(r) = p(r), (25) 

^0 

where A = /dx^ + jdy^ -|- jdz^ is the 3-dimensional Laplacian, and 

JF|A/(r)](k) = -|k|2 2^|/(r)](k) = -|k|q'(k), (26) 


where /(k) = J^[/(r)](k). As a result, the electrostatic potential of the point charge fl22|) has 
Coulomb’s form 


<h(r) 


1 Q 
47r£o |r| 


(27) 


The second well-known case [2211211 ES] is such that 


^||(|k|) — £o(^l + 


1 

r|,|kp 


(28) 
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Substituting fl28l) into (El]), we obtain 


( 29 ) 


(l'‘r + rr) «'k = ^Pk. 

Then, using the inverse Fourier transform of fl2TD . we get 

1 1 

A$(r)-9-*h(r) =-p(r). 

As a result, we have the screened potential of the point charge (j2^ in Debye’s form: 


$(r) 


1 Q 

AttSq |r| 


exp 





(30) 


(31) 


where is Debye’s radius of screening. It is easy to see that Debye’s potential differs from 
Coulomb’s potential by factor CD(|r|) = exp(—Irl/r^j). Such factor is a decay factor for 
Coulomb’s law, where the parameter ro dehnes the distance over which signihcant charge 
separation can occur. Therefore, Debye’s sphere is a region with Debye’s radius rm in which 
there is an influence of charges, and outside of which charges are screened. 


2 Fractional Power-Law of Non-Locality and General¬ 
ized Debye’s Screening 

In this section, we consider power-law type generalizations of Debye’s permittivity (l28ll . and 
generalizations of the correspondent equations for electrostatic potential <F(r) of the form 

A$(r)-w*h(r) = —^p(r) (32) 

by involving the fractional generalization of the Laplacian HD SUSIE]. 

2.1 A Power-law Generalizations 

In this section we consider the simplest power-law forms of the longitudinal permittivity e||(|k|) 
and correspondent equations for the electrostatic potential. The suggested simple models allows 
us to consider new possible types of an anomalous behavior of media with fractional power-law 
type of non-locality. We introduce some deformation of power-law type to well-known model 
of Debye’s screening. 

The generalized model is described by the deformation of two terms in equation fl28]) for 
permittivity in the from 

e||(|k|) = £o(|kr ^ (22) 
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The parameter a characterizes the deviation from Coulomb’s law due to non-local properties 
of the medium. The parameter jS characterizes the deviation from Debye’s screening due to 
non-integer power-law type of non-locality in the medium. 

Substituting fl3^ into we obtain 

(W" + T|k|dl>k = ^/>k. (34) 

and using the inverse Fourier transform of fl5T|) . we have 

((-A)“/=<I.)(r) + 4((-A)4/=*)(r) = i-p(r), (35) 

where (—and (—are the Riesz fractional Laplacian, see, for instance, 
and Appendix 1. Note that r and td are dimensionless variables. 

In order to describe the properties of two types of deviations separately, we consider the 
following special cases of the proposed model. 

1) Fractional model of non-local deformation of Coulomb’s law in the media with spatial dis¬ 
persion defined by equation (l33|) with (3 = 0, given by 

^||(|k|) = £o(|k|“"^ + 

Then equation fl?I|) has the form 

(|kr + 4)<^>k = ^Pk, (37) 

\ r^/ Sq 

and the equation for electrostatic potential is 

((-A)“-'=4.)(r) + T$(r) = lp(r). (38) 

This model allows us to describe a possible deviation from Coulomb’s law in the media with 
nonlocal properties defined by power-law type of spatial dispersion. 

2) Fractional model of non-local deformation of Debye’s screening in the media with spatial 
dispersion is defined by equation fl33D with a = 2, is given by 

^||(|k|) = £o(l + ij^p-p)' (39) 

Equation (1^ with (159]) lead to 

(|k|^ + i|k|'’)'S>k = )-ft. (40) 
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and then the corresponding equation for generalized potential is given by 


A$(r) + 4((-A)^^'<^>)(r) 


D 


1 


P(r), 


(41) 


which involve two different differential operators. Such model allows us to describe a possible 
deviation from Debye’s screening by non-local properties of the plasma-like media with the 
generalized power-law type of spatial dispersion. 

The behavior of electrostatic potentials for fractional differential models described by equa¬ 
tions fl38|) and flTT|) will be consider in Section 3.3. To the mentioned model can be hnd a 
explicit solution in terms of a Green type function. Also we will describe analytic solutions of 
the fractional differential equations fl5SD . 

2.2 General Fractional Power-Law Type of Non-Locality 

In the more general case, we can consider the following power-law form 


m 



(42) 


Substitution of (I42|) into and using the inverse Fourier transform gives the fractional 

partial differential equation 



(43) 


where am > ••• > cxi > 0, and aj G M (1 < j < m) are constants. 

We apply the Fourier method to solve fractional equation (1431) . which is based on the relation 


_F[(-A)”'y(r)|(k) = |kr/(k), 

Applying the Fourier transform to both sides of (1431) and using fl44)) . we have 


(44) 



(45) 


The fractional analog of the Green function (see Section 5.5.1. in [5]) is given by 
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where a = (ai, am)- 
The following relation 

r /'o_'jn/2 /-oo 

/ e-'‘‘’'>/(|k|)<i"k=-L^ / /(A)A"''2j„/2_,(A|r|)<iA (47) 


holds (see Lemma 25.1 of [1]) for any suitable function / such that the integral in the right- 
hand side of flT7|l is convergent. Here Jy is the Bessel function of the hrst kind. As a result, the 
Fourier transform of a radial function is also a radial function. 

On the other hand, using flT7|) . the Green function flT 6 l) can be represented (see Theorem 
5.22 in [5]) in the form of the one-dimensional integral involving the Bessel function J 1/2 of the 
first kind 


where we use n = 3 and a = (oi,..., am)- Note that 





+ flo I Ji/ 2 (A|r|) dA, 


(48) 


^ 1 / 2 ( 2 :) 



(49) 


If am > 1 and Am 7 ^ 0, Aq 7 ^ 0, then equation fl43|l (see, for example. Section 5.5.1. pages 
341-344 in [5]) has a particular solution is given by (|5n|l . Such particular solution is represented 
in the form of the convolution of the functions G(r) and p(r) as follow 

$(r) = — f G„(r — r') p(r') d^r', (50) 

^0 JR3 

where the Green function Ga(z) is given by fITSD . 

Therefore, we can consider the fractional partial differential equation (|43|) with oq = 0 and 
oi 7 ^ 0, when m G N, m > 1, and also the case where ai < 3, am > 1, > 1, cq 7 ^ 0, 7 ^ 0, 

am > ... > Q?! > 0 , which is given by 


Hi 

“T ^0 


(51) 


The above equation has the 


following particular solution (see Theorem 5.23 in |S]), given by 
$(r) = — f G„(r — r') p(r') d^r', (52) 

^0 JR3 


with 

|„|-l/2 roo / "i \ 1 

= (5(77? I ^’''"^i/2(A|r|)<lA, ( 53 ) 

These particular solutions allows us to describe electrostatic held in the plasma-like media with 
the spatial dispersion of power-law type. 


9 





2.3 Potentials for Particular Cases of Non-Integer Power-Law Type 
Non-Locality 


In this section we will study the properties of electrostatic potentials for fractional differential 
models mentioned in Section 3.1. 

Here we consider particular solutions of the following fractional partial differential equation 

((-A)"-'=4)(r) + a4(-A)'’'=4)(r) = ip(r), (54) 

where 1 < a, 0 < /5 < a, /9 < 3, and ap = Note that r and r^) are dimensionless. Equation 
fl5^ is the fractional partial differential equation fl5T]) with m = 1, and such equation has the 
following particular solution 

- r')p(r')c^^r', ( 55 ) 

^0 JR3 

where the Green type function is given by 

11-1/2 poo _ 

G„,s(r) = (|A|“ + as|A|«)-‘A=>/"J,/ 2 (A|r|)<iA. (56) 

Therefore, the electrostatic potential of the point charge 0221 ) has form: 


$(r) 


1 Q 
dvr^o li"! 


■c=p(kl). 


(57) 


with 


CaA\A) 


2 /-- Asin(A|r|) 
TT Jo |A|" + a^|A|/^ 


where Ga,^(|r|) describes the difference of Coulomb’s potential. 

Now we will study three special cases: 1) a 7 ^ 2, a > 1 and /9 
3) a 7 ^ 2 and (3 > 0. 


0; 2) a 


(58) 


2 and 0 < /5 < 2 ; 


2.3.1 Non-local deformation of Coulomb’s law (the case /5 = 0) 


Fractional model of non-local deformation of Coulomb’s law in the media with spatial dispersion 
is defined by equation (l5T|) with = 0, is given by 

((-A)“/2$)(r) + ^$(r) = —p(r), (59) 

A £0 

where a > 1, and oq = r~A > 0. Using 0571) and 0581) . it is easy to see that the electrostatic 
potential *h(r) differs from Coulomb’s potential by the factor 


Qo(|r|) 


2 r AsMAin) ^^ 

TT Jo |A|" ao 


(60) 


10 






Note that Debye’s potential differs from Coulomb’s potential by the exponential factor C£)(|r|) = 
exp(-|r|/rD). 


In Figure 1 we present plots of Debye exponential factor CD(|r|) = exp(—|r|/r£)) and gen¬ 
eralized factor Ca_o(|r|) for the different orders of 1.5 < a < 3.0 and oq = = 1. 


Using (Section 2.12.1 No.3. page 169. in |12]), we obtain the asymptotic (|r| —)■ 0) in the 
form 

22-“r((3-a)/2) 1 


C'«,o(|r|) ~ 
C'a,o(|r|) ~ 
Qo(|r|) 


\2-a 


-)2—a 


10-2 


^/^TT{a/2) I 
r((3-a)/2) 
■v/Fr(a/2) ' 

2r(3/a)r(l-3/a) 

1-3/q 
TTtt Oq 


, (1 < a < 2), 


, (2 < a < 3), 


r|, (a > 3). 


(61) 

(62) 

(63) 


Note that asymptotic fl6UI62p for 1 < a < 2 and 2 < a < 3 does not depend on the parameter 
Oq. We point out that for a = 2, using (Equation (11) of Section 1.2. in the book [B]), we 
obtain Debye’s exponents C 2 ,o(|i'|) = C'd(Ii'I). 

As a result, the electrostatic potential of the point charge in a media with this type of 
spatial dispersion will have the form 


$(r) 


Q 22-“r((3-a)/2) 1 


\/^r(a/ 2 ) 


13—0 


(1 < a < 2, 2 < a < 3) 


(64) 


on small distances |r| -C 1. In the case a > 3, we have the constant value of the potential for 
|r| <C 1 given by 

$(r) ^ T—(« > 3), (65) 

Aneo Reff 

where Reff is an effective sphere radius that is equal to 


-^e// ~ 


71 a a. 


1-3/q 

0 


2r(3/a)r(l-3/a)' 


( 66 ) 


Therefore, the electric field E is equal to zero at small distances |r| -C 1. It is well-known that 
the electric field inside a charged conducting sphere is zero, and that the potential remains 
constant at the value it reaches at the surface. Then the electric field of a point charge in 
the media with power-law of spatial dispersion with a > 3 is analogous to the field inside a 
conducting charged sphere of the radius Reff, for small distances |r| -C 1. 

The study of the asymptotic behavior of C'a^ 2 (|i'|) for |r| —)■ oo is an open question, although 
we have evidences to can suggest, as a conjecture, that its asymptotic behavior follow a power- 
law type also. Also from the corresponding plots, we can observe that the C'„^ 2 (|i’|) decreases 
more slowly than Debye’s exponent C£)(|r|). 
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It is easy to proof that C'q,^ 2 (|i’|) has a maximum for the case 2 < a < 3 and the maximum 
does not exists for 1 < a < 2, while for the particular case a = 2 it is well-known that it is the 
classical exponential Debye’s screening. 


2.3.2 Non-local deformation of Debye’s screening (the case a = 2) 


Fractional model of non-local deformation of Debye’s screening in the media with spatial dis¬ 
persion is described by equation fl5^ with a = 2, given by 

(67) 

has the following form 

( 68 ) 


- A$(r) + a^((-A)*'''^<S)(r) = —p(r), 

^0 

where 0 < /3 < 2. The electrostatic potential of the point charge 


■«>(■■) = 7^0 ■ c^A\A), 


where the function 




2 /■“ 

C2A\^\) = - 


A sin(A|r|) 


dA, 


(69) 


"■Jo + 

where C' 2 ,/ 3 (|r|) describes the difference of Coulomb’s potential 

In Figure 2 we present some plots of Debye exponential factor C£)(|r|) = exp(—|r|/r£)) and 
factor C 2 ,/ 3 (|r|) for different orders of 1.5 < /3 < 2 and ap = 1. 


Using (Equation (1) of Section 2.3 in the book 
behavior for C' 2 ,/ 3 (|r|) with /? < 2, when |r| —)• oo 


9 

C' 2 ,,(|r|) = - 


A sin(A|r|) 


TT Jq A^ -1- U/jA^ 


d\ ^ AqA)- 


, we obtain the following asymptotic 
1 ^ 1 


12-/3 




k=l 


.|(2-/3)(fc+l)’ 


where 


AoA) = — r(2 - /?) sin , 

Trap \2 / 


Am = 


Tra 


/c+l 

f 




(70) 

(71) 

(72) 


As a result, we have that generalized non-local properties deforms Debye’s screening such 
that the exponential decay is replaced by the following generalized power-law 

^0 




12-/3 


(0 < /3 < 2 ). 


(73) 


On the other hand, the electrostatic potential of the point charge in the media with this type 
of spatial dispersion is given by 

^0 Q 


$(r) 


A-KEn 


13-/3 


(0 < /? < 2 ) 


(74) 
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on the long distance |r| S> 1. 


2.3.3 Non-local deformation of Coulomb’s law and Debye’s screening (the case 
a ^ 2 and /9 > 0) together 


The electrostatic potential for non-local fractional differential model that is described by equa¬ 
tion includes two parameters {a,f3), where a > (3 > 0. In such model non-local properties 
deforms Coulomb’s law and Debye’s screening such that we have the following fractional power- 
law decay 


C'a,/3(|r|) ~ 


2r(2 -/3) sin(7r/j/2) 
nap 


1 

r 12-/3 


(|r| -» oo), 


(75) 


for 0 < /3 < 3 and a> (3. Note that this asymptotic behavior |r| —)■ oo does not depend on the 
parameter a. The field on the long distances is determined only by term with (—(a > (3) 
that can be interpreted as a non-local deformation of Debye’s (second) term in equation (l32ll . 

The new type of behavior of the spatial-dispersion media with power-law non-locality is 
presented by power-law decreasing of the field at long distances instead of exponential decay. 

The asymptotic behavior C'a^/jdrl) for |r| —)■ 0 is given by 


Q/^drl) ~ 

22-“r((3-a)/2) 1 

(1 < a < 2), 

(76) 

V^r(Q;/2) rp““’ 

CaA\^\) ~ 

2 =-“r(( 3 -a)/ 2 ) 

(2 < a < 3), 

(77) 

^/nT{a/2) 

CaAA) 

2 , , 

1—3/q; • fo ! \ ^ 

aa^ sin(d7r/aj 

(a > 3), 

(78) 

where we use Euler’s reflection formula for Gamma function. Note that the above asymptotic 
behavior does not depend on the parameter A ci-nd relations fl76tl77p does not depend on ap. 
The field on the short distances is determined only by term with (—A)“/^ (a > (3) that can be 


considered as a non-local deformation of Coulomb’s (first) term in equation fl32l) . 

On the other hand, it is remarkable that exist a maximum for the factor Co^/ 3 (|r|) in the 
case 0 < (3 < 2 < a. 


3 Fractional Weak Spatial Dispersion 

Here we introduce a generalization of well-known weak spatial dispersion for the power-law 
type of non-locality of the media [23l |2ll [25l [26l |23 1^ - 
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3.1 Weak Spatial Dispersion 

Let us give an short description of weak spatial dispersion in the plasma-like media (for details 
see, for instance, [23l [2ll |25] ). 

Spatial dispersion in electrodynamics is called to the dependence of the tensor of the absolute 
permittivity of the medium on the wave vector [221 IMl ES] • It is well-known that this depen¬ 
dence leads to a number of phenomena, for example the rotation of the plane of polarization, 
anisotropy of cubic crystals and other [26l [271 [28l |29l [30l [3TI [32l [33l [HI [35l |36l l37] . 

The spatial dispersion is caused by non-local connection between the electric induction D 
and the electric field E. Vector D at any point r of the medium is not uniquely defined by the 
values of E at this point. It also depends on the values of E at neighboring points r', located 
near the point r. 

Non-local connection between D and E can be understood on the basis of qualitative anal¬ 
ysis of a simple model of the crystal. In this model the particles of the crystal lattice (atoms, 
molecules, ions) oscillate about their equilibrium positions and interact with each other. The 
equations of oscillations of the crystal lattice particles with the local (nearest-neighbor) inter¬ 
action gives the partial differential equation of integer orders in the continuous limits [33 SO]- 
Note that non-local (long-range) interactions in the crystal lattice in the continuous limit can 
give a fractional partial differential equations [39l 00]. It was shown in [39l 00] that the equa¬ 
tions of oscillations of crystal lattice with long-range interaction are mapped into the continuum 
equation with the Riesz fractional derivative. 

The electric field of the light wave moves charges from their equilibrium positions at a 
given point r, which causes an additional shift of the charges in neighboring and more distant 
points r' in some neighborhood. Therefore, the polarization of the medium, and hence the held 
D depend on the values of the electric helds E not only in a selected point, but also in its 
neighborhood. This applies not only to the crystals, but also to isotropic media consisting of 
asymmetric molecules and plasma-like media [221122] [22] . 

The size of the area in which the kernel djj(r) of integral equation fiT2]) is signihcantly 
determined by the characteristic lengths of interaction i?o- For different media these lengths 
can vary widely. The size of the area of the mutual inhuence Ro are usually on the order of 
the lattice constant or the size of the molecules (for dielectric media). Wavelength of light 
A is several orders larger than the size of this region, so for a region of size Rq value of the 
electromagnetic field of light wave does not change. By other words, in the dielectric media 
for optical wavelength A usually holds kRo ~ i?o/A ~ 10“^ <C 1. In such media the spatial 
dispersion is weak [231 Ell EHl EZ] • To analyze it is enough to know the dependence of the tensor 
£jj(k) only for small values k and we can replace the function by the Taylor polynomial 


^u(k) ^ij T 'yijiki -|- ^ijimkik^ -|- ... . (79) 

Here we neglect the frequency dispersion, and so the tensors 'jiji, 5ijim do not depend on 
the frequency u. 

The tensors in fl72]) are simplihed for crystals with high symmetry [22]. For an isotropic 
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linear medium, we can use 


e||(|k|) = £ + 7|k| + 5|k|^ + ... . (80) 

In order to explain the natural optical activity (for example, optical rotation, gyrotropy) is 
sufficient to consider the linear dependence on k in fl7^ and fl80l) . For non-gyrotropic crystals 
it is necessary to take into account the terms quadratic in k. 

For power-like type of non-locality we should use fractional generalizations of the Taylor 
formula (see Appendix 2). 


3.2 Fractional Taylor series approach 

The weak spatial dispersion in the media with power-law type of non-locality cannot be de¬ 
scribes by the usual Taylor approximation. The fractional Taylor series is very useful for 
approximating non-integer power-law functions |38]. To illustrate this point, we consider the 
non-linear power-law function 

e||(|k|) = a„|k|"-F oq. (81) 

If we use the usual Taylor series for the function flHTj) then we have inhnite series. 

For fractional Taylor formula of Caputo type (see Appendix 2), we need the following known 
property of the fractional Caputo derivative (see, for instance, 0) 

aDk{k-af = (x > a, a > 0, f3 > 0), (82) 

i (p — a -I- 1) 

where k = |k|. In particular, if /3 = a, then 

-a)» = r(a+l), = 0. (83) 

Therefore 

(^B”£||)( 0 ) = r(a + 1 ), 

while, the higher order Caputo fractional derivatives of £||(|k|), given in (IHTD . are all zero. 
Hence, the fractional Taylor series approximation of such function is exact. Note that the order 
of non-linearity of £||(|k|) is equal to the order of the Taylor series approximation. 


3.3 Weak spatial dispersion of power-law types 


We consider such properties of the media with weak spatial dispersion that is described by the 
non-integer power-law type of functions £||(|k|). The fractional differential model is used to 
describe a new possible type of behavior of complex media with power-law non-locality. 

The weak spatial dispersion (and the permittivity) will be called cx-type, if the function 
£||(|k|) satishes the condition 


lim 

|kK0 


g||(|k|) -g||(0) 

|k|“ 


CLfy , 


(84) 
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where a > 0 and 0 < |aQ,| < oo. Here the constant Eq is the vacunm permittivity (^o ~ 
8.854 10-12^-m-i). 

The weak spatial dispersion (the permittivity) will be called (a,/3)-type, if the fnnction 
£||(|k|) satisfies the conditions fl8T)) and 


e:||(|k|) - £||(0) - aa£o |k|“ 
l™o £o|k|/5 


0/35 


(85) 


where (3 > a > 0 and 0 < |a/ 3 | < cxo. 

Note that these definitions are similar to dehnitions of non-local alpha-interactions between 
particles of crystal lattice (see Section 8.6 in [TS] and [321 SO]) that give continnons medinm 
eqnations with fractional derivatives with respect to coordinates. 

For the weak spatial dispersion of the (a,/9)-type, the permittivity can be represented in 
the form 

^||(|k|) = £o(^ + Oo|k|“ + a^|k|^) -|- i?Q:^^(|k|), (86) 

where £ = e:||(0)/£o can be considered as the relative permittivity of material, and 


lim 

|kK0 


|k|/5 


0 . 


As a result, we can use the following approximation for weak spatial dispersion 


(87) 


ir|l(|k|)/eo M £ + ao|k|“ + 


( 88 ) 


If a = 1 and /3 = 2, we can use the usual Taylor formula. In this case we have the well- 
known case of the weak spatial dispersion [26l [271 12H1 1291 [301 |3T1 |32l [33l [331 |35l [Ml EZ] In 
general, we should use a fractional generalization of the Taylor series (see Appendix 2). If the 
orders of the fractional Taylor series approximation will be correlated with the type of weak 
spatial dispersion, then the fractional Taylor series approximation of e:||(|k|) will be exact. In 
the general case, /3 ^ a, where 0 < /3 — a < 1, we can use the fractional Taylor formula in the 
Dzherbashyan-Nersesian form (see Appendix 2). For the special cases /3 = 2a, where a < 1 
and/or (] = a + 1, we could use other kind of the fractional Taylor formulas. In the fractional 
cases new types of physical effects may exist. 


3.4 Fractional differential equation for electrostatic potential 

We can consider a weak spatial dispersion of the power-law (a, /3)-type. Then substituting (j88il 
into CD, we obtain 

(e\k\‘^ + aa|k|“+^ = —Pk, (89) 

where £ = £||(0)£o and (3 > a > 0. The inverse Fourier transform of fl89l) gives 

a/ 3 ((-A)(^+ 2 )/ 2 ^)(^^ ^ a„((-A)(“+2)/2^)(r) _ ^^^(r) = -p(r). (90) 

^0 


16 




This fractional differential eqnation describes a weak spatial dispersion of the (a,/5)-type. 
Eqnation fITO has the following particular solution 

$(r) = — [ G2,aA^ - !■') pA) (91) 

^0 Jr3 

where G 2 ,q ,/3 is the Green function of the form 

11 - 1/2 foo 

G2,a,/3(r) = ^27r)3/2 J I1I ^ I1^) ^ Ji/2 (A| r |) dA. (92) 

Therefore, the electrostatic potential of the point charge fl22|) for this case is given by 

KW^ T^n-Cw(|r|), (93) 

4vr£o |r| 

where 0 < a < (3, and the function 


C2,aAA) = 


A sin(A|r|) 


d\ 


TT Jo aa|A|"+2 + a;3|A|^+2 + £|A|2 
describes the difference between such generalized potential and Coulomb’s potential. 
For the weak spatial dispersion of the a-type, we have 


(94) 


a„((-A)(“+=>''"*)(r) - £A<*(r) = -p(r). 

^0 


(95) 


This equation is a special case of equation fl90|) . where = 0. Then, the electrostatic potential 
of the point charge has form 

1 n 

(96) 


= ■C2,a,„(|r|), 


dvr^n |r| 


where 0 < a < A and 


C2,a,o(|r|) - - 


A sin(A|r|) 


d\. 


(97) 


TT Jo a„|A|“+2 + £|AP 

This case is described by the fractional differential model introduced in Section 3 for the 
case of the order of Riesz fractional derivative is a + 2 > 2. 

To describe properties of electric field of the point charge in the media with weak spatial 
dispersion, we consider properties of the function 


C2,aAA) = 


A sin(A|r|) 


Jo a/ 3 |A|^+^ + aQ,|A|“+^ + £|A|' 


dX, 


(98) 


where j3 > a. 
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Using the values for the sine integral Si{x) for the infinite limit 


sinu 


71 

2 ’ 


(99) 


and the equation for the integral transform (Section 2.3, equation (1) in [H]) of the form 


z°‘ ^ sin( 2 ;) = r(a) sin 


71(3 


, (-1 < a < 1 ), 


( 100 ) 


we obtain the asymptotic for the function 6 * 2 , 0 ,/?(II'D of the form 




/7ra\ i 


r(/3)sm(^^j). (101) 


It allows us to obtain the asymptotic behavior of the electrostatic potential 

Q 1 UaQ 2r(Q;) sin(7rQ;/2) 1 2r(/5) sin(7r/9/2) 1 


$(r 


47reeo |r| 47re^ec 


7T 


|a+l 


477e‘^en 


TT 




+ 1 ' 


( 102 ) 


where Q < a < (3 < 1. 

The parameter e is interpreted as a relative permittivity of the media. It is well known that 
far from the electric dipole the electrostatic potential of its electric field decreases with distance 
|r|, as |r|“^ (see Section 40 in [15]), that is faster than the point charge potential (*h(r) ~ |r|~^)- 
The first term in (110211 describes the well-known Coulomb’s field. The second and third 
terms in (I102p look like changed dipole electrostatic field that for integer case (a = 1) has the 
from <h(r) = dcos 0 /(dvr^o |rp), where d = |d|, and d is the (vector) dipole moment, and 9 is 
an angle between the vectors d and r (see Section 40 in [15]). We can consider the effective 
values 

7 f \ _ 2r(Q))flQ,Q _ ^/I \ I'l 

for non-integer values of a. The second and third terms in equation ( 11021 ) can be interpreted 
as a generalized dipole fields of power-law type with the non-integer orders a and (3, and these 
terms are represented as 


^e//(r) 


deffja) cos 9eff{a) _ deffi/3) cos OeffiP) 
dvr^o |r|"+^ dvr^o |r|^+i 


(104) 


where 0 < a < (3 < 1. 

In Figure 4 present some plots (see Figure 3) of Coulomb’s electrostatic potentials ‘h(|r|) = 
(l/dvr^o)and the potential with factor Cp^a, 2 {\'^\) for different orders of 0 < a < /? < 2 , 
where = e = 1 . 


From the plots (see Figure 4) it is easy to see that far from the point charge in the media 
with weak spatial dispersion the electrostatic potential decreases with distance |r| more slowly 
(11021) than the potential of the point charge potential (*h(r) ~ |r|~^)- 
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4 Conclusion 


We consider fractional power-law type generalizations of permittivity, and generalizations of 
the correspondent eqnations for electrostatic potential $(r) by involving the fractional gen¬ 
eralization of the Laplacian HB 021 n ig. The simplest power-law forms of the longitndinal 
permittivity £||(|k|) = are snggested. The parameter a characterizes 

the deviation from Conlomb’s law dne to non-local properties of the medinm. The parameter 
(3 characterizes the deviation from Debye’s screening dne to non-integer power-law type of non- 
locality in the medinm. The correspondent eqnation (13^) for electrostatic potential <h(r) that 
has the form ((—A)“/^<h)(r) -|-r^^((—A)^/^<h)(r) = £Q^p(r), contains (—A)“/^ and (—A)^/^ are 
the Riesz fractional Laplacian 5110210115], and r and td are dimensionless variables. 

To the mentioned model can be hnd a explicit solntion in terms of a Green type fnnc- 
tion. Also we will describe analytic solntions of the fractional differential eqnations 0351) 
for electrostatic potentials. The electrostatic potential of the point charge has form $(r) = 
Q/(47r£o|r|) • Ga^/jdrl), where Ga,/3(|r|) is dehned by fl58l) describes the differences of Conlomb’s 
potential and Debye’s screening. Using the analytic solntions of the fractional differential 
eqnations for electrostatic potentials, we describe the asymptotic behaviors of the electrostatic 
potential. The new type of behavior of the spatial-dispersion media with power-law non-locality 
is presented by power-law decreasing of the held at long distances instead of exponential decay. 

In order to describe the properties of deviations separately, we consider the following special 
cases of the proposed model: 

1) Fractional model of non-local deformation of Coulomb’s law in the media with spatial 
dispersion that corresponds to the case /? = 0 and a 7^ 2, a > 1. This model allows us to 
describe a possible deviation from Coulomb’s law in the media with nonlocal properties dehned 
by power-law type of spatial dispersion. 

The electrostatic potential of the point charge in a media with this type of spatial dispersion 
has the form <F(r) ~ |r|"“^ for 1 < a < 2 and 2 < a < 3 on small distances |r| 0. In the 

case a > 3, we have the constant value of the potential for |r| —)■ 0. Therefore the electric held 
E of a point charge in the media with power-law type of spatial dispersion with a > 3 is equal 
to zero at small distances |r| —)• 0 that is analogous to the well-known case of the held inside a 
conducting charged sphere of the radius Re//, for small distances. The asymptotic behavior of 
potential for |r| —)■ 00 follow a power-law type also by our assumption. From the corresponding 
plots, we observe that the Ca,2(|i'|) decreases more slowly than Debye’s exponent Ci^drl). The 
function 6*0^2(|r|) has a maximum for the case 2 < a < 3 and the maximum does not exists 
for 1 < a < 2, while for the particular case a = 2 it is well-known that it is the classical 
exponential Debye’s screening. 

2) Fractional model of non-local deformation of Debye’s screening in the media with spatial 
dispersion is dehned by a = 2 and 0 < /? < 2. Such model allows us to describe a possible 
deviation from Debye’s screening by non-local properties of the plasma-like media with the 
generalized power-law type of spatial dispersion. 

The generalized non-local properties deforms Debye’s screening such that the exponential 
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decay is replaced by the fractional power-law, and the electrostatic potential of the point charge 
in the media with this type of spatial dispersion is given by <h(r) |r|^ ^ for 0 < /3 < 2 on 

the long distance |r| —)■ cxo. 

3) Fractional non-local model that is described by equation fl5^ includes two parameters 
{a, f3), where a > f3 > 0 such that a 2. In such model non-local properties deforms Coulomb’s 
law and Debye’s screening such that we have the fractional power-law decay for 0 < /3 < 3 and 
a > (3. The asymptotic behavior |r| oo does not depend on the parameter a. The held 
on the long distances is determined only by term with (—that can be interpreted as a 
non-local deformation of Debye’s term. It is remarkable that exist a maximum for the factor 
in the case 0 < < 2 < a. 

The asymptotic behavior for |r| —)■ 0 is given by $(r) ~ |r|"“^ for 1 < a < 2 and 2 < a < 3, 
and $(r) is a constant for a > 3. Note that the above asymptotic behavior does not depend 
on the parameters jS and r]^. The held on the short distances is determined only by term with 
(—that can be considered as a non-local deformation of Coulomb’s term. 

We also consider weak spatial dispersion in the media with fractional power-law type of non¬ 
locality. In general, it cannot be describes by the usual Taylor approximation. The fractional 
Taylor series is very useful for approximating non-integer power-law functions. The media 
with spatial dispersion is described by the non-integer power-law type of functions £||(|k|). 

Using fractional generalization of the Taylor series (see Appendix 2), we get approximations 
of the form e:||(|k|)/£o ~ £ + aa|k|" + a^|k|^ for such type of media. If a = 1 and (3 = 2, 
we have the usual Taylor formula, and the well-known case of the weak spatial dispersion 
[261 \27[ [28l |29l [301 [SH |32l [33l |3ll [35l |36l [37] . In general, we should use a fractional generalization 
of the Taylor series. If the orders of the fractional Taylor series approximation will be correlated 
with the type of weak spatial dispersion, then the fractional Taylor series approximation for 
£||(|k|) will be exact. 

These fractional weak spatial dispersions is described by equation of the form a/3((—A)(^+^^A<|))(];-)_|_ 
aa((—— eA<h(r) = —p(r). To this fractional differential equation we can be 
hnd a explicit solution in terms of a Green type function. Also we describe analytic solu¬ 
tions of the fractional differential equations for electrostatic potentials. It allows us to obtain 
the asymptotic behavior of the electrostatic potential of the form <F(r) ^Couiomb{r) + 

^^^\dipoiei^) + ^e%,dipoiei^)^ where 0 < a < (3 < 1. The hrst term describes Coulomb’s 
held. The second and third terms are interpreted as electrostatic helds generalized of changed 
dipoles of power-law type with the non-integer orders a and (3, which have the form $g“y(r) = 
—deff{a) cos 9eff {a)/ where 0 < a < (3 < 1, deff is the effective dipole mo¬ 
ment, and 6 is an effective angle. For integer case {a = 1) we have the usual from <F(r) = 
dcos9/{ATreo |rp) of the helds. 
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Appendix 1: Riesz fractional derivatives and integrals 

Let us consider Riesz fractional derivatives and fractional integrals. The operations of frac¬ 
tional integration and fractional differentiation in the n-dimensional Euclidean space M" can 
be considered as fractional powers of the Laplace operator. For a > 0 and ’’sufficiently good” 
functions f{x), x G M”, the Riesz fractional differentiation is dehned [H] 021 0115] in terms of 
the Fourier transform by 


(-A);/V(;r) = ;F-'(|k|“(.F/)(k)). (105) 

The Riesz fractional integration is dehned by 

n/(i) = 7^-‘(ik|-“(7^/)(k)). (106) 

The Riesz fractional integration can be realized in the form of the Riesz potential 011021 0115] 
dehned as the Fourier convolution of the form 

= / Ka{x - z)f{z)dz, (a > 0), (107) 

where the function Ka{x) is the Riesz kernel. If a > 0, and a ^ n,n + 2,n + A, the function 
Ka{x) is dehned by 

If a 7 ^ n, n -f 2, n -7 4,..., then 


K^{x) = --fj{a)\x\°‘ "^ln\x 
The constant 7n(tt) has the form 


7n(a) 


2"7r”/^r(Q(/2)/r(^^) aj^n + 2j, n G N, 

(_l)(n-«)/22«-i7rV2 r(a/ 2 ) r(l + [a - n]/2) a = n + 2j. 


(108) 


Obviously, the Fourier transform of the Riesz fractional integration is given by 

7-(n/(i)) = ikr“(7^/)(k). 

This formula is true for functions f{x) belonging to Lizorkin’s space. The Lizorkin spaces of 
test functions on is a linear space of all complex-valued inhnitely diherentiable functions 
f{x) whose derivatives vanish at the origin: 


T = {/(x) : /(x) G ^(M"), {D^fm = 0, |n| G N}, (109) 
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where S'(M”) is the Schwartz test-function space. The Lizorkin space is invariant with respect 
to the Riesz fractional integration. Moreover, if f{x) belongs to the Lizorkin space, then 

rj{x)iif{x) = n+'’/n), 


where a > 0, and /3 > 0. 

For a > 0, the Riesz fractional derivative (—= —d°‘/d\x\°' can be dehned in the form 
of the hypersingular integral (Sec. 26 in [1]) by 




1 

dn{m,a) 


1 

2-\ct-\-n 


{ATf){z)dz, 


where m > a, and {A^f){z) is a finite difference of order m of a function f{x) with a vector 
step z G M"' and centered at the point x G M"': 


{ATmz) = Y,(-^y- 

j=0 


m\ 


j\ {m-j)\ 


f{x 


jz). 


The constant dnim, a) is defined by 


_7r^+-/M^(a)_ 

2"r(l -|- Q!/2)r(n/2 -|- a/2) sin(7rQ!/2) ’ 


where 




Note that the hypersingular integral (—A)“'^^/(x) does not depend on the choice of m > a. 

If f{x) belongs to the space of ’’sufficiently good” functions, then the Fourier transform 
of the Riesz fractional derivative is given by 


(7-(-A)“'V)(k) = |k|“(7^/)(k). 


This equation is valid for the Lizorkin space |1] and the space of inhnitely differentiable 

functions on M"" with compact support. 

The Riesz fractional derivative yields an operator inverse to the Riesz fractional integration 
for a special space of functions. The formula 

(-Ar/‘rj[x) = !(xX (a> 0 ) ( 110 ) 


holds for ’’sufficiently good” functions /(x). In particular, equation flllUp for /(x) belonging 
to the Lizorkin space. Moreover, this property is also valid for the Riesz fractional integration 
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in the frame of Lp-spaces: f{x) G Lp{M.) fori ^ p < n/a. Here the Riesz fractional derivative 
(—is understood to be conditionally convergent in the sense that 

= lim(-A);/^ (111) 

c_:».n ’ 


where the limit is taken in the norm of the space Lp(M), and the operator (—A)x!e 
by 


(-A) 


a/2 

X,€ 


1 

dn{m,a) 



1 


{ATf){z)dz, 


is dehned 


where m > a, and {A^f){z) is a hnite difference of order m of a function f{x) with a vector 
step 2 ; G M"' and centered at the point x G M”. As a result, the following property holds. If 
0 < a < n and f{x) G Lp{R) fori ^ p < n/a, then 


(-ArJH“f(x) = f{x), (OO), 


where (—is understood in the sense of (I1 11 1I . with the limit being taken in the norm of 

the space Lp(M). This result is proved in [1] (see Theorem 26.3). 

We note that the Riesz fractional derivatives appear in the continuous limit of lattice models 
with long-range interactions |18] . 


Appendix 2: Fractional Taylor Formula 

Riemann-Liouville and Caputo derivatives 

The left-sided Riemann-Liouville derivatives of order a > 0 are dehned by 

1 f d\^ r f{x')dx' 


RL D(A 




T{n — a) \dxJ {x — 
We can rewrite this relation in the form 

dx 


, (n = [al + 1). 


(“B„V)(x)= T- (4”r/)(!). 


( 112 ) 


(113) 


where /", is a left-sided Riemann-Liouville integral of order a > 0 


(4“+/)(^) = 


fix') dx' , , 

7 TT —1 > a). 

/H—a ’ ^ ^ 


r(«) A ix-x') 

The Caputo fractional derivative of order a is dehned by 

d 


(114) 


GBsV)M = (4V ( ^ I /IW. 


(115) 


26 










where is a left-sided Riemann-Liouville integral (llldh of order a > 0. In equation (11271) we 
use 0 < a < 1 and n = 1. The main distinguishing feature of the Caputo fractional derivative 
is that, like the integer order derivative, the Caputo fractional derivative of a constant is zero. 

Note also that the third term in fll27p involves the fractional derivative of the fractional 
derivative, which is not the same as the 2a fractional derivative. In general. 

Then the coefficients of the fractional Taylor series can be found in the usual way, by repeated 
differentiation. This is to ensure that the fractional derivative of order a of the function (x —a)“ 
is a constant. The repeated the fractional derivative of order a gives zero. Then the coefficients 
of the fractional Taylor series can be found in the usual way, by repeated differentiation. 

Fractional Taylor series in the Riemann-Liouville form 

Let /(x) be a real-valued function such that the derivative (^'^Zl"+™'/)(x) is integrable. Then 
the following analog of Taylor formula holds (see Chapter 1. Section 2.6 i): 

/(^) = Y1 rf ‘"1 ■ ~ + Rm{x), (a>0), (116) 

T{a + 3 + l) 

where are left-sided Riemann-Liouville derivatives, and 

R„.{x) = ( 117 ) 

Riemann formal version of the generalized Taylor series 

The Riemann formal version of the generalized Taylor series [ISIIIT]: 

-hOD 

/(^) = 

m=—c 

where for a > 0 is the Riemann-Liouville fractional derivative, and = Ia°‘ for 

a < 0 is the Riemann-Liouville fractional integral of order |a|. 

Fractional Taylor series in the Trujillo-Rivero-Bonilla form 

The Trujillo-Rivero-Bonilla form of the generalized Taylor formula [15] : 

m 

fix) = Y{{j + l)a) + Rmix, a), (119) 


RLpa+mf^^Xo) 
r(a -I- m -I- 1) 


(x - Xo) 


a-\-m 


(118) 
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where a G [0; 1], and 


( 120 ) 

( 121 ) 


Cj = r(Q;) |(x - a)' “('* D“Yf(x)](a+), 

Rm(.x,a) J) (a: “) . 

Fractional Taylor series in the Dzherbashyan-Nersesian form 

Let ttfc, {k = 0,1, be increasing seqnence of real numbers such that 

0<afc-afc-i<l, ao = 0, /c = 1, 2,m. (122) 

We introduce the notation [l9l |50] (see also Section 2.8 in 0): 

(123) 


^0+ 


In general, Dq^. Fractional derivative differs from the Riemann-Liouville 

derivative by finite sum of power functions since (see Eq. 2.68 in |^) 


ra 7-/3 / T 

-'o+-'o+ W -'0+ 


The generalized Taylor formula 


(124) 


m—1 


f{x) = ^akX°‘>^+Rm{x), {x>0). 


k=0 


where 


(D^ 

r(a» + l) ’ ' 


r(Q;m + 1) Jo 




(125) 


(126) 


Fractional Taylor series in the Odibat-Shawagfeh form 

The fractional Taylor series is a generalization of the Taylor series for fractional derivatives, 
where a is the fractional order of differentiation, 0 < a < 1. The fractional Taylor series with 
Caputo derivatives [51] has the form 


f{x) = f{a) + 


r(Q! + 1) 


(x - a)“ + 




T{2a + 1) 


(x - a)^" + ..., 


(127) 


where is the Caputo fractional derivative of order a. 
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Figure 1: Plots of Debye exponential factor Cd[x) = exp(—with 
= Ca,o{x) with ao = = 1 for the orders: a) a = 1.9, b) o; = 1.6, c) 

ere x = Irl and we use 0 < a; < 5. 








Figure 2: Plots of Debye exponential factor Cd{x) = exp(—with = 1 and the factors 
y = C 2 ^ii{x) with ap = 1 for the orders: a) /9 = 1.9, b) /9 = 1.6, c) /9 = 1.1, d) /9 = 0.6. Here 
X = Irl and we use 0 < a; < 5. 
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Figure 4: Plots of Coulomb’s electrostatic potential <F(a;) = x~^ and the electrostatic potentials 
y = x~^ ■ C 2 ,a,i 3 {x) (the weak spatial dispersion case) with = e = 1 for the orders: a) 

a = 0.1 and ^ = 1.1, b)Q; = 0.2 and {3 = 0.4, c) a = 0.4 and (3 = 1.4, d) a = 0.8 and (3 = 1.8. 
Here x = |r| and we use 0 < a; < 100. The asymptotic behavior for long distances is dehned 
by 13 only. 









